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Abstract. 

In this article we discuss the Dirac equation in the presence of an attractive 
cylindrical <5-shcll potential V(p) = —aS(p — po), where p is the radial coordinate and 
a > 0. We present a detailed discussion on the boundary conditions the wave function 
has to satisfy when crossing the support of the potential, proceeding then to explore 
the dependence of the ground state on the parameter a, analyzing the occurrence of 
supercritical effects. We also apply the Foldy-Wouthuysen transformation, discussing 
the non-relativistic limit of this problem. 
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1. Solution of the Dirac equation for a cylindrical 5-shell 

The Dirac equation opened certainly a new world to physics, with marvelous and 
surprising consequences as, for example, the existence of anti-matter, the fact that 
g = 2 for the electron giromagnetig moment, etc. [l]-[3]. However, perhaps one of the 
most interesting facts associated with the Dirac equation is the idea of the existence of 
a populated vacuum, the Dirac sea, i.e. a vacuum state having a non trivial structure. 
This concept developed later into a cornerstone of our present understanding of modern 
quantum field theory. Remnants of this exciting idea are, for example, the vacuum 
condensates in Quantum Chromodynamics, the thermal vacuum distribution, when 
finite temperature effects are taking into account in heavy ion-collisions, and, in general, 
the idea of non-zero vacuum expectation values. Supercritical effects, i.e. an instability 
that appears when the ground state starts to dive into the Dirac sea, inducing then 
positron emission, is a prediction of this scenario, which has nevertheless not been 
experimentally confirmed yet [I]. 

Singular 5-type potentials have been considered in quantum mechanics already from 
the beginning. In general, singular 5-type potentials can be considered as toy models 
which allow us to get a physical insight, being, at the same time, more easy to deal 
with in comparison with more realistic extended potentials. An example is the well 
known Dirac-comb or Kronig-Penney model [5] in non relativistic quantum mechanics, 
which provided us with an understanding and intuition about the emergence of band 
structures in solid state physics. Another situation where this type of interactions have 
been used, is in nuclear physics as a model for residual interactions between nucleons 
inside an incomplete nuclear shell, [6], providing a correct qualitative picture to this 
difficult problem. 

In the frame of the Dirac equation, singular 5-type interactions have also been 
considered many times in the literature [7], [8]. Certainly this problem is also attractive 
from the perspective of mathematical physics. A rigorous construction of self adjoint 
extensions for the Dirac operator, allowing the handling of matching conditions at 
the support of the 5-potential for the spherically symmetric case were considered in 
[9] . These conditions are essential for determining the existence of bound states. The 
theorem by Svendsen (TO] is also a rigorous result which tells us that for the spherically 
symmetric singular potential V(p) = —a5(p — p ), supercritical effects will be absent 
in the limit where po — ► 0. This fact is a particular case of the general result of this 
theorem, namely that contact interactions for the Dirac equation can be constructed only 
if the supporting manifold has a codimension of dimension one, which is the case in our 
problem. See also [TT] in this context. The cylindrical attractive V(p) = —a 5(p — po) 
potential, the main subject of this article, as we will see, presents many interesting 
features, being worthwhile to go into a detailed discussion of its properties. This case 
has been discussed in the literature but in the context of the (2 + 1) Dirac equation, 
where we have three two by two Dirac matrices [12J. The treatment here is different, 
and relies on an appropriate unitary transform and the use of the chiral representation 
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for the Dirac matrices. 

If we think about the interaction of a neutral Dirac particle, which carries a 
magnetic moment /i, for example a neutron, with an external magnetic field B, 
Hj = —ft ■ B, we will have a point-like interaction of a 5- type if the external magnetic 
field corresponds to an Aharonov-Bohm magnetic vortex. The situation we propose to 
discuss in this article corresponds physically to an array of Aharonov-Bohm vortices, 
distributed along a circle of radius po, which may interact with the magnetic moment 
of a neutral Dirac particle. Several experiments have been done on the behavior of 
electrons under the influence of external Aharonov-Bohm magnetic fields [13], and such 
an array could in principle be constructed. 

Obviously, we are compelled to work in cylindrical coordinates. Since a 5-type potential, 
which has support in a domain of zero measure, divides the space into two regions, we 
will start by considering first the free particle case. Then, in a second step we will 
establish the connection between the wave functions when crossing the support of the 
potential. 



1.1. Dirac free particle in cylindrical coordinates 

The Dirac equation in curvilinear coordinates, in general, includes a non trivial spin 
connection in the covariant derivative [H]. When considering the free particle case 
in cylindrical coordinates, it turns out that the relevant Dirac matrices are coordinate 
dependent because here we have 

( e~ i9 \ ( -ie~ il 

"> ~ \f ) 5 ae ~ [ie i6 

This is something we would like to avoid when solving the free Dirac equation 
(7 • p — m)ip = 0. To deal with this, let us consider the unitary transformation 

S=±eW, (2) 
VP 

introduced in [15]. Applying this transformation to the Dirac equation and defining 

ip 1 = Sip, we get 

( 7 °<9 t + 7^ + jd e + 7 3 <9 Z + imW = 0. (3) 

Notice that in our case, due to the form of 5", the spin connection term is trivial. At this 
point we can see the effect of the unitary transformation S. It has rotated our equation 
so that the 7 matrices involved are the cartesian ones, eliminating their dependence on 
the coordinates. 

We now introduce the Ansatz ip' = 7 1 7 2 0, and rewrite equation (3) as 




(if! + £ 2 ) = 0, 



(4) 



M. Loewe, F. Marquez and R. Zamora 



4 



where 

Hi = (7% + 7 3 ^) 7 V (5) 

H 2 = (V^+^ + w^V- (6) 

It is easy to see that these operators commute. This means that if each one of them 
satifies an eigenvalue problem, they have a common basis, i.e. 

H x <j> = \ l( j> (7) 
#20 = A 2 0. (8) 

Putting this into equation (4), it is trivial to note that Ai = — A 2 = A. In this way we 
can write two equations for the 4-component bispinor </> 

(-7 2 <9 P + -jd e + im^i 2 + \)<j> = (9) 

(lVl 2 d t + 7 3 7 V<9, - A)0 = 0. (10) 
It is natural then to propose 

= e -iEt e ik z z e ik e m ^ ^ 

where e and rj are both 2-component spinors. Note that, because of the symmetry of 
the problem we would expect to have a plane wave solution in the z coordinate. This 
ansatz is quite general, precisely, due to the strong symmetry of our problem. Since we 
are imposing no restriction on the radial coordinate p, this ansatz should represent a 
general and complete solution of the problem. For our purpose it is convenient to use 
the Chiral representation for the Dirac matrices [16] given by 





/= „ , y= " . (12) 



Substituting then our Ansatz in the previous equations we get 
iEo^o 1 ^ + ik z a i o 1 o 2 ri - Xe = 

iEo^o^e - %k z o i (j 1 o 2 e - Xrj = 0. (13) 

It's important to emphasize that even though we have a rotated spinor and we 
are working in the Chiral basis, our energy levels remain unchanged since unitary 
transformations do not affect the energy spectrum. We define 



v= r , d4) 




where now, e lj2 and rji^ are scalar complex functions (i.e. each of them represents one 
of the components of the bispinor 0). Replacing this on the previous equation, we can 
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obtain the following relations 

A 2 = E 2 - k\ 
E + K 

£1= ^ — m (is) 

E - k z 

These relations are quite simple. Their derivation would have been more involved if 
the normal Dirac representation for the Dirac matrices would have been used. From 
equations (9) and (10), using the previous relations, we get 



k e .(m + X)(E + k Z) 



V2 = V2 + i : Vi (16) 

P A 

/ k B , (jn - \){E - k z ) 

Vi = —m - * r V2, (17) 

P A 

where a prime denotes derivative with respect to p. 

Without loss of generality, we can take k z = decoupling this set of equations. We 
define r = £p, where £ 2 = m 2 — A 2 . It is important to notice that ^ 2 is positive definite, 
since we are looking for bound states, i.e. — m < E < m. In terms of this new variable 
we can write 



r 2 V2 ~ {r 2 + k e (k e + = 0, (18) 

The previous equation is the well-known Bessel equation of second kind. In this way, 771 
and t]2 are given by 

7] 2 = p^(AI ke+l/2 (£p) + BK ke+1/2 (£p)) 

1 

m = 4 C^j) 2 ph-AI ke -i/2(£p) + BK ke . 1/2 {ip)) (19) 

= ixpH-Ahe-idtp) + BK ke _ l/2 (ip)). (20) 

Once we have found the solution for the Dirac free particle in cylindrical coordinates, 
we can return to our original problem, the attractive 5-shell potential. The boundary 
conditions when crossing the support of the 5 potential, to be discussed in the next 
section, will provide the energy spectrum. 

1.2. Boundary conditions for the Dirac equation in the presence of a cylindrical 5 -shell 

By applying the transformation S to the Dirac equation, (j-p — m + ^°a5(p — po))ip = 0, 
we get 

( 7 °ft + 7^ + ^d e + 7 3 «9 Z + im - ia^°5(p - p W = 0. (21) 
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kg 

<^> - id p T] 2 - i — 772 - imei - a5(p - p )??i + Aei = (22) 
P 

kg 

id p rji - i—Vi + i m£ 2 + a $(p - Po)V2 + Ae 2 = 0. (23) 

From here we get 

1 d f n 2 \ . r/ s 
— +iab{p - po) 



rfi/rfi -I dp \rji 



i—ViV2 + im(5i?7i - e 2 r] 2 ) - K £ iVi + £ 2V2) ) • (24) 



v 2 - vi V P 



Integrating between the limits [p — e , Po + e ] ; we have 



( ~ ) 



To keep going, we divide the space in two regions: 
Region I: p < p 



Vx = -ixp^AI k0 „ l/2 (26) 
r? 2 = p^AI kg _ 1/2 . (27) 



Region II: p > p 



Vi = -ixP^BK kg _ 1/2 (28) 
ri 1 = p*BK ke _ 1/2 . (29) 

The functions in regions I and II are different because the Dirac four-spinor has to be 
normalized. In this way we get the relation 

vl ' (po+e) vi(p °- €) ---tana. (30) 



r i2 I (po+ e ) \ ( vi(po-t) 



1 



We are interested in the ground state energy (kg = 0). Using the explicit expressions of 
the relevant Bessel functions [18], and taking the limit e — > 0, we get the transcendental 
equation that determines the behavior of the ground state 

i=f + tanh (s ) 
i+ e .-t. ~ n = - tan a > ( 31 ) 



1- tanh (s ) 



where we have defined the following dimensionless variables 
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E 
m 

ip = mp (32) 
s = (fiVl - e. 

For the derivation of this equation it was important to fix the value of the constant A. 
Its value can either be ±|i£|. In order to get rid of this ambiguity we notice that when 
a = (no potential), the energy of the electron should correspond to the free particle 
energy. To accomplish this A has to be equal to — \E\. Finally, we can remove any 
ambiguity by demanding a continuous energy spectrum when E — 0. This demands 
A = —E. 

The evolution of the ground state energy as function of the different parameters involved 
is shown in Fig.l. Numerically it seems that there will be no critical effects, i.e || — > 
when e — > — 1. In fact it is easy to show this property analytically. If we consider eq. 
(31) and define 1 + e = e, then we need to consider the limit when s — > 0. It is easy to 
find for the leading behavior of the derivative 



da cos 2 a(tanh(A/2mp ) + tan a) 
when e — > 0, i.e. supercritical effects in this case are absent, in contrast to what occurs 
in the case of the spherical delta potential discussed in [7]. 




Figure 1. The energy of the ground state for a <5-shell potential as a function of the 
coupling a. 
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2. Foldy-Wouthuysen transformation for the <5-shell 

We would like now to discuss the Foldy-Wouthuysen [T7j transformation in the case of 
our 5-shell problem. In this section we will obtain some transcendental equations for 
the energy of the ground state up to order 1/m and 1/m 2 separately. Later, we will 
compare this energy with the energy provided by the Schrodinger equation. In addition, 
we will show how the spin-orbit coupling term appears from the Foldy-Wouthuysen 
transformation up to order — 



2.1. Foldy-Wouthuysen transformation for the Hamiltonian 
We have the hamiltonian 

H = a ■ p + f3m — a5(p — po). (33) 

We apply a first Foldy-Wouthuysen trasformation F = e lS , with S = —ip^, and then 
a second one with S' = — (£o ~ 2m a ^(P ~ Po)) where, similarly to the previous one, 
all odd operators in the hamiltonian have been included inside S'. We also have defined 

-*2 

£o = 111 this manner we arrive to 



H" = (3^m+^j -aS(p-p ) 



a 



]m 2 



(2(5 -p6(p- Po))a g p e - a ■ p (a ■ p 5(p - p ))) (34) 



We can see that H" has no odd operators, and so it should be the non relativistic limit of 
the Hamiltonian up to order 1/m 2 . It is not difficult to see that, in spite of the presence 
of two a matrices, the last two terms in equation (34) correspond to even operators. 

2.2. Solving the free particle in the non relativistic limit 

Now we solve the free particle using the hamiltonian H" without the Dirac delta 
potential. We start from 

M =0 ( (35) 



dt \ 2m r 
and suggest the ansatz ip = e lEt \ I with u = e lke u(p), to get 



v 



r 2 u" + ru' - k 2 e u - k 2 u = 0, (36) 

where n 2 = —2m(E — m) and r = up. 

It is important to note that, since we are looking for bound states, i.e. — m < E < 
m, k 2 is positive definite. 
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Equation (36) is the Bessel equation, and its solutions are well known. We can 
follow a completely similar process for v to get the solutions 

u = AI kg (np) + BK kg (np) (37) 
v = CI kg ( ^2m{E + m)p) + DK kg ( ^2m{E + m)p) . (38) 

Now that we have the solutions, we can apply the boundary conditions, associated to 
the S potential. 

2.3. Boundary conditions for the non relativistic limit of the 5 -shell 
We have the Dirac equation, up to order 1/ra 

'lit = if m + ^ ~ a5{p ~ Po) ) ^ (39) 
For the upper component, we get 

— k 2 u = u" H — v! — -^u + 2ma5(p — po)u. (40) 
P P 

We are interested in the ground state (kg = 0), so the previous equation becomes 

— k 2 u = u" H — u + 2ma5(p — po)u. (41) 
P 

We integrate between [po _ e, Po + e ]- We should keep in mind that, since the equation is 
a second order differential equation, unlike the usual Dirac equation, our wave function 
is continuous in p : 

u'Qt e e + 2mau(p ) = 0. (42) 

Once again, we divide the space in two regions. 
Region I (p < p ): 

Ul = AI q {k P ) (43) 
Vl = CIo ( v / 2m(E + m)p) . (44) 
Region II (p > p ): 

u n = BKo(Kp) (45) 
v n = CK (^2m(E + m)p). (46) 

These functions are different in each region in order to assure normalization of the Dirac 
spinor. 

With this, we can write, taking the limit e — > 

u'jj -u'j + 2mau(p ) = 0. (47) 
Imposing continuity of u in p we can get the relation 

A = B K f p0 }. (48) 
h(Kp ) 



M. Loewe, F. Marquez and R. Zamora 10 

Using this result and defining the dimensionless variables 

_ E 
m 

ip = p m. (49) 

we can write 



V^l-^i^poKo^Po) 



^/2(1 -e)K (Kp )h(K Po ) - 2aK (K Po )I (K Po ) = 0, (50) 



which is a transcendental equation for the energy of the ground state, up to order 1/m. 
We now repeat the previous procedure up to order 1/m 2 . We have the Dirac equation 

%d t = (z 3 ( m + L) - a5{p - po) 

(2(3 -p5{p- po))a e p e - a ■ p (a ■ p 6{p - p ))) ) ^, (51) 



8m 2 
where 

2(d-p5(p- p ))a e pe = -25' (p - p )a p a e W e = -26' (p - Po )a z V e . (52) 

If we focus in the upper component, this term becomes: 

- 25\p - Po )cr z Ve = -A5\p - Po )L z S z . (53) 

We recognize the spin-orbit coupling, which is one of the relativistic corrections provided 
by the Foldy-Wouthuysen transformation. However this term has no influence on the 
ground state, where L z = kg = 0. We now proceed to solve equation (51). Using the 
same ansatz as before for the upper component, we get 

7-1 XL u kn , . 

Eu = mu- h -u - ab{p - p )u 

2m 2mp 2mp A 

' k > a > y(fi - Pa)u -^( y y- p „ ) + lk^A) u . (54) 



4m 2 p 8m 2 \ p 

We solve for the ground state (kg = 0) and the previous equation becomes 
pu" v! 



2m 2m 



pa6(p - p )u 



om \ p 



Integrating between [p — e, p + e] and taking the limit e — > 0, we get 



7T-( u 'ii(Po) ~ u'i(po)) + au(po) + -^< 7 (p ) + ^2«/j(a>) = 0. (56) 
2m Sm z sm z 



M. Loewe, F. Marquez and R. Zamora 11 
After inserting solutions (43) and (45), equation (56) becomes 

- a/2(1 - £)I q (kpo)K 1 (kpo) + (a + ~(1 - e)J I (kpo)K (kpo) 

-)J^p-Ii(Kpo)K o ( K po) = 0. (57) 

This is yet another transcendental equation for the energy of the ground state, up to 
order 1/m 2 . 

From equations (57) and (50), we can see how the ground state energy depends on 
the coupling constant a. 




Figure 2. Behaviour of the ground state as a function of a. The solid line is the 
solution to Schrddinger's equation and the dashed line is the second order Foldy- 
Wouthuysen approximation, both for tp = 1. 

Although we have not mentioned it here, the solid line stands for both the 1/m 
order of the non relativistic limit and the Schrodinger equation, (after the rest energy has 
been considered), so there is no relativistic correction, other than the rest energy, to the 
Schrodinger equation at first order in the Foldy-Wouthuysen transformation. However, 
up to order 1/m 2 , the first relativistic corrections begin to appear. As we discussed, the 
spin-orbit coupling can be explicitly written in the equation. Other corrections must be 
present to, since the energy of the groundstate decreases for the same coupling value, 
with respect to the non relativisitic solution. 

3. Conclusions 

We have discussed the solutions of the Dirac equation for an atractive cylindrical delta 
potential. It turns out that there are no supercritical effects, in contrast to what happens 
in the spherical symmetric case. We discussed also the non relativistic limit of these 
problems according to the Foldy-Wouthuysen approximation. No relativistic corrections 
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are found up to order 1/m, however, up to order 1/m 2 relativistic corrections are present, 
such as the spin-orbit coupling. 

4. Acknowledgements 

The authors acknowledge support from FONDECYT under grant Nr. 1095217. M.L. 
acknowledges also support from the Proyecto Anillos ACT119. F.M. and R.Z. 
acknowledge support from CONICYT under grants 21110577 and 21110295. The 
authors would like to thank M. Banados for helpful discussions about the spin 
connection. 

References 

[1] P.A.M. Dirac, Proc. Roy. Soc. (London) A117, 610 (1928). 
[2] P.A.M. Dirac, Proc. Roy. Soc. (London) A118, 351 (1928). 

[3] B. Thaller, The Dirac Equation, Texts and Monographs in Physics, 1992 (Berlin: Springer). 

[4] See W. Greiner, B. Miiller, and J. Rafelski, Quantum Electrodynamics of Strong Fields, Text and 

Monographs in Physics, 1985 (Berlin: Springer); B. Miiller, Ann. Rev. Nucl. Sci. 26 351 (1976); 

W. Greiner and J. Reinhardt, Quantum Electrodynamics, ch. 7, 2009 (Berlin: Springer); W. 

Greiner and J. Reinhardt, Theoretical aspects of quantum electrodynamics in strong fields, Vol. 

440, p. 153, 1994 (Berlin: Springer). 
[5] R. de L. Kronig and W. G. Penney, Proc. Roy. Soc. (London) A 130 (1931) 499. 
[6] See, for example, J. D. Walccka, Theoretical and Subnuclear Physics, second edition, World 

Scientific and Imperial College Press 2004. 
[7] R. D. Benguria, H. Castillo, and M. Loewe, Journal of Physics A: Mathematical and General 33, 

5315 (2000); M. Loewe and M. Sanhueza, Journal of Physics A: Mathematical and General 23, 

553 (1990). 

[8] D. Villarrocl, European Journal of Physics 19, 1998, 85. 

[9] J. Dittrich, P. Exner, and P. Scba, Journal of Mathematical Physics 30 (1989) 2875. 
[10] E. C. Svendsen, Journal of Mathematical Analysis and Applications 80 (1981) 551. 
[11] F. A. Coutinho and Y. Nogami, Physical Review A 42, 5716 (1990). 
[12] Shi-Hai Dong and Zhong-Qi Ma, Foundations of Physics Letters 15, 171 (2002). 
[13] See, for example, M. Peshkin and A. Tonomura, The Aharonov-Bohm Effect, Lecture Notes in 

Physics 340, 1989 (Berlin: Springer). 
[14] M. Nakahara, Geometry, Topology and Physics, fourth edition, section 7.10.2, Institute of Physics 

Publishing, Bristol and Philadephia, 1995; P. Schliiter, K.H. Wietschorke and W. Greiner, J. 

Phys. A: Math. Gen. 16 1999 (1983). 
[15] V. M. Villalba, Nuovo Cimento 112 B 109 (1997). 
[16] M.N. Hounkonnou and J.E.B. Mendy, J. Math. Phys 40 4240 (1999). 
[17] L.L. Foldy and S.A. Wouthuysen, Phys. Rev. 78, 29 (1950). 

[18] M. Abramowitz and LA. Stegun, Handbook of mathemathical functions: with Formulas, Graphs, 
and Mathematical Tables, New York: Dover (1965). 



